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Abstract of the Dissertation 


Concentration Susceptibility and the Phase Diagram of 
He 3 - He 4 Mixtures near the Tricriticai Line 

by 

Frederick Hermann Wirth 
Doctor of Philosophy 

in 

Physics 

State University of New York at Stony Brook 
June, 1978 

The variation of He 3 concentration with height is measured in 
the lower phase of He 3 -He 4 mixtures near the tricriticai point 
(TCP), at pressures of 2.5, 9.9, and 16.7 atmospheres. A 
moveable plunger is used to vary the height of the phase boundary 
relative to a tunnel diode oscillator concentration sensor. The 
sample is isolated at low temperature to maintain concentration 
integrity. Relative shifts of 0.0002 in concentration and 0.06 mm. 
in height are discernable. 

The variation is shown to fit a functional form of Ax e = Z where 
x is the relative change in concentration from its critical value 
and Z is the vertical distance below the interface. When x is 
expressed in mole% and Z in centimeters, the parameters take 
on the values given below at the pressure investigated. 
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cm* 


Pres sure (atm.) 



16.7 ±0.1 


2.5 ± 0.1 


9.9 ±0.1 


4.7 ± 0.9 


6.8 ±2 


7.3 ±5 


1.94 ± 0.2 
1.81 ±0.48 
2.10±0.2 


The thermodynamics of He 3 -He 4 mixtures and the scaling 
hypothesis are discussed. The relationship between the 
gravitationally induced concentration gradient and the concen¬ 
tration susceptibility, (aX/SA) TP , is derived. It is shown that 
6 is the critical exponent defined by (X-JY C ) 6 = (A- A c ). Values 
of (9X/&A) Tp are calculated from the data and compared with 
values measured at saturated vapor pressure by other re¬ 
searchers. 

The theoretical models of van Leeuwen and Cohen, Blume, 
Emery and Griffiths and Takagi are discussed. The values 
they predict for 6 are presented. An attempt is made to under¬ 
stand the pressure dependence of the concentration susceptibility 
in the light of these models. Improved determinations of the 
scaling properties of these models near the TCP are 
briefly discussed. 

Values of the concentration gradient are used to provide 
gravitational corrections to data of Johnson and Rohde. Cor¬ 
rected data showing the motion of the TCP between saturated 
vapor pressure and 22 atm (the tricritical line) are presented. 
The determination in this work that the tricritical concentration 






















takes on its maximum value at 9.9 ± 0,4 atm is incorporated in 
this presentation. Comparisons are made between the pressure 
dependence of: the T = 0K solubility of He 3 in He 4 , the value of 
X c , and the magnitude of the concentration susceptibility. 
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I. Introduction 


Isotopic mixtures of He 3 and He 4 display some very unusual 
behavior in the liquid phase. The experimentally accessible 
thermodynamic variables for the investigation of this behavior 
are: the temperature - T; the pressure -P; and the molar con¬ 
centration of He 3 -X. Upon cooling, provided that X *£ 0.67, the 
entire mixture will become superfluid at some temperature , T x . 
This higher-order, superfluid phase transition is known as the 
“lambda” transition because of the characteristic shape of the 
associated specific heat curves. T x is a monotonically de¬ 
creasing function of X which maps out a line in T-X space 
known as the “lambda line”. 

Upon further cooling, provided that the mixture 

undergoes a first-order phase transition, separating into two 
liquid phases of different X. The upper phase, rich in He 3 floats 
on the denser, He 4 -rich phase. The He 4 -rich phase is super- 
fluid . This phase transition can also be displayed in T - X space 
as a "coexistence curve”, which displays, for each value of T, 
the coexisting concentrations appropriate to the upper and lower 
phases. 

Both the lambda line and the coexistence curve typical of a 
He 3 -He 4 mixture at saturated vapor pressure (SVP) are dis¬ 
played in Fig. 1. In order to demonstrate the use of this diagram, 
consider a mixture of average concentration, X = 0.5. Following 












down the solid line marked ( p ), as temperature is 
lowered, we first encounter a lambda transition at 7\ = 1.3 K 
at which the entire sample becomes superfluid. Upon further 
cooling, the sample begins to phase separate at T s -0.8 K into 
phases of concentrations X,. and X 2 . At some lower temperature, 
T = 0,4 K, phases of concentrations X s and X A coexist in relative 
amounts given by the ratio of the lengths L 2 and L t respectively. 
Note that at T-OK some He 3 remains dissolved in the He 4 . This 
is a result of the Fermi statistics of He 3 and is the basis for the 
operation of the He 3 dilution refrigerator. 

Many systems are known that display one or the other of these 
two types of phase transition. Higher-order transitions dis¬ 
playing similar specific heat curves occur in binary alloys. 1 
Many mixtures of organic liquids, including phenol and water, 
have a range of temperatures over which they are not mutually 
soluble in all proportions. However, systems that display both 
phenomena in this manner are quite rare. In particular, only a 
few anisotropic antiferromagnetic materials such as FeC^ and 
dysprosium aluminum garnet 3 have phase diagrams similar to 
those of He 3 -He 4 mixtures. 

It is reasonable to infer from Fig. 1 that these two phase tran¬ 
sitions are, in fact, closely related phenomena. The coexistence 
curve shows an angular peak atX = 0.67 at which point the lambda 
line terminates exactly. This intersection divides the coexistence 
curve into the so called left-hand and right-hand branches. The 












Figure 1 


The He 3 - He 4 phase diagram under SVP 
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usual coexistence curves for binary mixtures have rounded tops 
much closer to-?f = 0.50. 4 It will come as no surprise, then, 
when phenomenological model calculations for this system indi¬ 
cate that the same interactions which cause the superfluid 
transition are responsible for the phase separation as well, 

A region of particular interest on this phase diagram is the 
point of intersection of the lambda line and the coexistence curve, 
known as the tricritical point (TCP). At this point, both the 
phase transitions described above occur simultaneously and con¬ 
tinuously . All of the data in this work were taken in the im¬ 
mediate neighborhood of this point, in order to elucidate some 
of the rather complex phenomena that occur there. 

We will later show that the quantity (dX/d&) TP diverges at 
this point (A = /i 3 - p 4 , the difference in chemical potential be¬ 
tween He 3 and He 4 ). This quantity is analagous to a magnetic 
susceptibility in that it is a measure of how sensitive the con¬ 
centration of the mixture is to any change in its driving “field”. 
Therefore the sample becomes infinitely "soft” at the TCP and 
any change in A which is a function of T t P , and.X', will cause 
infinite shifts in concentration. For this reason, a close ap¬ 
proach to the TCP is a matter of some delicacy. The major 
contributor to inhomogeneities in concentration as one approaches 
the TCP is gravity. Previous workers 5 ’ 13 have estimated concen¬ 
tration gradients in excess of 1% over a height of 1 cm. in the 
lower phase at the TCP. Since the exact magnitude and the form 
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of this gradient were not known, it became difficult to extend the 
left and right hand branches of the coexistence curve into the im¬ 
mediate neighborhood of the TCP. Indeed, the concentration at 
a given height in the sample chamber might follow a path in T-X 
space like the one shown by the line of points in Fig. 1, skirting 
entirely around the TCP and not joining the true coexistence 
curve for 10 mK below the TCP. 

The original motivation for this work was to determine the 
magnitude and the functional form of the gravitationally induced 
concentration gradient near the TCP in order to determine the 
phase diagram more precisely and completely. We simultaneous¬ 
ly hoped that this information would also make possible a 
determination of (9X/9 A) tp to sufficient precision to make a 
quantitative statement about its limiting behavior as the TCP is 
approached. 

This asymptotic behavior can be described in the scaling hy¬ 
pothesis (to be discussed later) by a parameter 5 defined as; 



( 1 . 1 ) 


Here we introduce the subscript c as indicative of the value of 
a variable at the TCP. 6 is the “critical exponent" associated 
with the singular behavior of the concentration susceptibility 
near this point. Some justification for the case of Eq. (1.1), as 
well as proof of the relationship between the gravitationally 
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induced concentration gradient and (DX/3 A) tp will be given in 
section m. 

The remainder of this work will be organized as follows: In 
section II we will consider the prior experimental investigation 
of the region near the TCP. In section HI a number of phenom¬ 
enological models will be presented. Their predictions will be 
compared with the previously discussed experimental results, 
and efforts will be made to extract the predicted trieritical be¬ 
havior for (5X/&A) T7 . In section IV the apparatus employed in 
this work will be described. The procedures used to extract data 
as well as the data itself, will be presented and analysed in section 
V- Section VI is reserved for a summary and conclusion. 










II. Previous Experimental Work 

The complex behavior of He 3 - He 4 mixtures has been under in¬ 
vestigation for nearly thirty years, Abraham, Weinstock and 
Osborne 7 observed the lambda line fromX^O to X = 0.28 in 1949 
by measuring the flow rate of a sample through a superleak as a 
function of temperature. The existence of a coexistence curve 
was demonstrated by Walters and Fairbank 8 in 1956. They em¬ 
ployed a nuclear magnetic resonance technique to observe the 
He 3 line with sensing coils at the top and bottom of their appara¬ 
tus. Subsequently, both the lambda line and the coexistence curve 
were explored by a variety of methods. Daunt and Heer 41 used the 
heat leak through a sample filled capillary to locate the superfluid 
transition while Dash and Taylor 42 employed an oscillating disk 
viscosimeter for the same purpose. Peshkov and Zinov’eva 43 
observed the coexistence curve optically while Roberts and 
Sydoriak obtained data from observation of discontinuities in the 
temperature derivative of first sound velocities. De Bruyn 
Ouboter et al. 45 determined the coexistence curve by using 
specific heat measurements. While this list of experiments is 
certainly not complete, it gives some idea of the wide variety of 
methods available to the experimenter in this field. The shape of 
the top of the coexistence curve and its intersection with 
the lambda line was first accurately determined by Graf, 

Lee and Reppy 3 in 1967. They used a capacitative 
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technique identical to the one employed in this experiment to 
determine the concentrations of the upper and lower phases. 

The lambda transition was detected by qualitative measurements 
of the thermal relaxation time. Many other experimenters have 
contributed greatly to this field of knowledge. A review of this 
early work can be found in Taconis and DeBruyn Ouboter 10 and 
Peshkov 

Subsequently, several experiments were done which mea¬ 
sured the effects of pressure upon the phase diagram and 
explored the region near the TCP in great detail. Notably 
Zimmermann 12 et al . employed specific heat measurements at 
SVP to determine Cp* , C p A, A and {BX/d A) tp over the ranges 
0.53 0.73 and 0.6 1.4 K. Goellner and Meyer 13 mea¬ 

sured the vapor pressure of He 3 -He 4 mixtures over a similar 
range of temperature and concentration. From this data they 
obtained values for A and (8A/8 A) tp and investigated the critical 
behavior of these quantities near the TCP. Watts and Webb 0 ' 14 
used light scattering from the critical opalescence of He 3 -He 4 
mixtures near the TCP to directly determine {aj£'/SA) TP . Again, 
the critical behavior of this quantity was investigated. The 
gravitationally induced concentration gradient also was mea¬ 
sured by an interferometric technique, and was determined to 
have resulted in a decrease in A' of 0.005 ±0.002 at a distance 
4mm below the interface between the two phases at T e , Abler s 
and Greywall 15 ’ 16 have used careful second sound measurements 












in the superfluid region near the TCP and along the left hand 
branch of the coexistence curve. The phase diagrams deter- 
mined by these researches all agree tolerably well. Both 
Goellner and Meyer, and Zimmermann ei at, also determine that 
the exponent 6 has a value close to 2. 

Away from SVP, Lai, 17 Johnson 18 and Hohde 19 have explored the 
phase diagram over the range 0.35 =X =0.80 at pressure varying 
from SVP to 22 atm.-near the solidification point-using the 
same techniques as Graf, Lee and Reppy, A composite of their 
data is shown in Fig. 2. Note that the essential feature of the 
phase diagram remains; the lambda line stays firmly rooted in 
the angular peak of the coexistence curve. The TCP itself shows 
interesting behavior, T c decreasing monotonically with increasing 
P while X c first increases with increasing P up to ~9.9 atm. and 
then decreases with increasing P until the solidification point is 
reached. It should be noted that the T = 0 solubility of He 3 in He 4 
as measured by Watson, Richardson and Reppy 30 shows qualita¬ 
tively similar behavior. 

Near the TCP the results of all these experiments show the 
effect of the presence of gravitationally induced concentration 
gradients. Attempts have been made to correct for the 
'‘broadening” of thermodynamic functions and the deviations in 
concentration by estimating - or in the case of Watts and Webb, 
directly measuring-the gradient. 















Figure 2 


The pressure dependence of the He 3 -He 4 phase diagram 
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m. Theoretical Treatments of He 3 -He 4 Mixtures 


Three different theoretical approaches to understanding the 
behavior of He 3 -He 4 mixtures near the TCP will toe considered 
here. They are: thermodynamic relations, the scaling hypoth¬ 
esis and critical exponents and finally, the predictions of 
several phenomenological models of the phase transition. Be¬ 
tween them, we shall be able to predict the form of the 
gravitational gradient, relate it to important thermodynamic 
quantities and gain some understanding of the microscopic pro¬ 
cesses that combine to produce it. 

A. Thermodynamic relations 


In order to include the concentration of our binary mixture in 
a thermodynamic formulation, let us begin with the usual form of 
the first law: 

dU - TdS - PdV + ^ 3^3 + ^dN A (3.1) 

where \i ( and N { are respectively the molar chemical potential 
and the number of moles of He*. We then define X as the molar 
concentration of He 3 in the mixture: 




N s +N 4 


If N=Ns+N 4 then 


(3.2) 


12 









(3.3) 


N^=NX 
N € =N(1~X) 

So that: 

dU - TdS - PdV + (p 3 - H 4 )NdX (3.4) 

If we divide through by N , and define a chemical potential 
difference by: 

A = (h 3 “^ 4 ) (3.5) 

we arrive at: 

du - Tds - Pdv + A dX (3. 6 ) 

where lower case variables indicate quantities per mole. 

Note the close analogy with magnetic systems in which a term 
HdM appears in the first law. In this case we have a “field" A 
that drives changes in concentration rather than a field H doing 
the same for the magnetic moments of a system. 

Let us now define the molar Gibbs potential of the system: 

G^u-Ts+Pv (3.7) 

So that: 

dG=-sdT+vdP + MX (3.8) 

This Gibbs potential must be a concave up function of concen¬ 
tration at all physically accessible points i.e.: 











(3.9) 


because a negative value for {B^/Bx)^ leads to instabilities in 
concentration just as a positive value of compressibility makes 
mechanical equilibrium impossible. We can similarly show that 
the Helmholtz free energy is concave up, i.e. that (aA/a*) TV ^0. 

Since we know that a phase separation does occur in He 3 -He 4 
mixtures, we may infer that the Gibbs potential looks something 
like Fig. 3, For T>T ct G(x) must be concave up for all values 
ofX. For T<T C there must be a "forbidden region” corre¬ 
sponding to the interior of the coexistence curve, for which 
G(x) is concave down. This is the region between B and C on the 
lower curve. The concentration of the upper and lower phases is 
found by a bitangent construction (solid line) which minimizes 
the total Gibbs energy of the system. From the nature of this 
construction it is apparent that: 



(3,10) 


So A is the same in both coexisting phases. Further, as T ap¬ 
proaches T c from below, the hump in the G(x) curve must get 
smaller and smaller, at last becoming a point of inflexion when 
X l ~X a , that is, at the TCP. This implies that: 



(3.11) 














Figure 3 


A possible Gibbs potential for He 3 - He 4 mixtures 




























Thus the concentration susceptibility, which is the inverse of the 
above quantity, diverges at the TCP. 

In order to include gravity in our system, we must take a 
closer look at the p 3 and p 4 in Eq. (3,1). Applying a method due 
to Griffiths , 21 let us define p^(T,P,^) as the molar chemical 
potential of He 1 free of the influence of gravity. Then; 

ti i (T,P t X,z) = n ( l{T t P J X) + m t gz (3.12) 

where m i is the gram molecular weight of He*, g is the accelera¬ 
tion of gravity and z is the vertical height from some arbitrary 
reference point. Then: 

A(T, P,X,z) -A° + {m 3 - mjgz (3.13) 

with: 

We can now begin to understand the mechanism that causes 
the appearance of the gravitational gradient. In equilibrium, A 
must be a constant throughout the sample. Were it not, it would 
force concentration changes in the sample, thus changing ^{TPX) 
until it became constant. Since A is constant in eq. (3.13) we have 
to assume that A 0 varies with height in such a way as to make the 
right hand side of that equation constant as well. As A 0 has no ex¬ 
plicit z dependence, its variation must be due to changes in X 
and P; T being a constant throughout the mixture at equilib¬ 
rium. Thus including gravitational energies in our formulation 










of the first law has forced the appearance of a gravitational 
gradient in concentration. 

Expressing the above idea in mathematical form: 

_/8A\ /0A\ dP fd A\ dX ,„ iA . 

dz U* /tpx + ^ p )xxz dz + W) TPZ dz (3 ' 5 

Tliis equation gives us a relationship between the concentra¬ 
tion susceptibility (a^/SA^ and dX/dz, the gravitational 
gradient, since we can evaluate or estimate ail other terms in 
the expression: 

From eq. (3.13): 


(: 


i7j =(m 3 -w 4 )^ 

V* / T px 


(»)„-( 


&^\ 

5Xj 


(3.15) 


= where p = the density of the mixture* (3.16) 
dz 

and from a Maxwell relation (Appendix B): 


(3.17) 


TP 


giving: 


( 


dX\ dX/dz _ 

8iJ TP g[ (m 4 ~mj) + p{dv/BX% p ] 


(3.18) 


Note that we have supressed the "z-held-constant” subscript 


^We will show in Appendix A that the approximation of constant 
density is a very good one. 







for clarity in this and subsequent thermodynamic relations. It 
will be understood that all subsequent thermodynamic deriva¬ 
tives will be taken at constant z or, equivalently, have their 
gravity-free values. The denominator of the right hand side of 
eq. (3.18) contains only well behaved functions. Therefore the 
critical behavior of the concentration susceptibility is mirrored 
in the behavior of the gravitationally induced concentration 
gradient. 


B. The Scaling Hypothesis 


The initial justification for the scaling hypothesis is experi¬ 
mental, It was noticed that sufficiently close to a critical point, 
the pairwise relationships between many suitably reduced ex¬ 
perimental quantities could be expressed as power laws. For 
instance, if we choose the reduced temperature, / 



(3.19) 


as our independent variable, and express some other thermo¬ 
dynamic quantity as /(/), then the limit: 



(3.20) 


gives us the critical exponent X . The implication is that suffi¬ 
ciently close to T c 




(3.21) 










This would be discovered experimentally by measuring / as 
a function of t close to t = 0 and plotting the result on a log-log 
graph. If the data lay on a straight line as t — 0, then the slope 
of this line would be the critical exponent for those two variables. 
An example of this process is shown in Fig. 4 in which the re¬ 
lationship between reduced density and temperature is explored 
near the gas-liquid critical point of He 4 . 

In more formal mathematical terms, the scaling hypothesis 
can give us the functional forms of (9 X/&A) tp and, hence, 
dX/dz as we near the TCP. We first define a generalized 
homogeneous function: 

F(jc } y) is homogeneous if: 


F(X a x,X b y) = XF(^,y) 


(3.22) 


F(x,y) is said to be homogeneous in x to degree “a” and 
homogeneous iny to degree u b’\ All homogeneous functions may 
be transformed by a suitable choice of X to: 


F(x,y) = x°f'{y/x b ); X = i 

A 


(3.23) 



(3.24) 


These relations will be useful in defining critical exponents. 
The scaling hypothesis then becomes the assumption that near 
the point of interest, the singular parts of the appropriate ther¬ 
modynamic potential are generalized homogeneous functions of 








Figure 4 


Measurements of the coexistence curve for the liquid-vapor 
phase transition in He" near its critical point. After Roach 38 
(1968), The critical exponent £ = 0.354 
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the thermodynamic fields. In particular, let us assume that the 
potential: 


A=u~T a+ Pv-AX (3.25) 

is homogeneous. (Then all other potentials related to this one by 
Legendre transformations are also homogeneous.) From eq, (3.6) 
and eq. (3.25) 


dA~- sdT 4 vdP - XdA 
Going over to the reduced variables: 
t = T-T e ; x^X-Xj A'=A-A c 
We have 


(3.26} 


(3.27) 


dA~— sdt 4 vdP - xdA ' 

LetrfP = 0. Then, by our assumption 
G(X a /,X*A')=XC(l,A') 


and: 


&G(\*t,\*A r ) _ 9G(f,A') 
SA' &A' 


(3.28) 


(3.29) 


(3.30) 


so: 


\>x{\ a t,\ 1 >A')=\x{t,A f ) 


(3.31) 


if we set X = /" l/c and explore the limit as / —0 along i'=0 we get 















{t’ L t*y^ x (l,0)=x{t,Q) 


(3.32} 


and we see that: 

lim x(f, 0 )~ f ( 3 , 33 ) 

t—Q 

A **0 

by the usual definition of critical exponents 

&~X c )cc |r- T c | fl (3.34) 

and we can make the identification £={1 - b)/a. 

Similarly, by setting X = A r ' L ^ b ) and proceeding as above to 
determine the behavior as A' — 0 along f = 0 we can get the rela¬ 
tion; 

limjf(f,0)-A ' 1 (3.35) 

i *-*0 

(iO 

and comparing this to the ■usual: 

(X'-Jf c ) <t cc(A-A c ) (3-36) 

we see: 

6 = —^— (3.37) 

1-b 

By choosing other values for X and taking other differentia¬ 
tions, the other usual critical exponents may be expressed in 
terms of a and 5, and thus, finally, in terms of one another. The 
results of these calculations, as well as the definitions of the 





usual critical exponents are summarized below in table {3,1), 

Table 3.1 

Definition 

©„'l r - r »l' T 

(x-^,)oc|T-r c | fl 
(A-A c )cc pf~X c > 
or 

cp - T (H) TO “i r -^i-“ 

y=J3(5-l) 
a+zfi+y-Z 

Another justification for the scaling hypothesis is that many 
of the relationships between the critical exponents, such as those 
in table 3,1, above, can be derived as inequalities by rigorous 
thermodynamic argument. Also, as we shall see, many theories 
give us thermodynamic results that "scale". 

We are now in a position to say something about the functional 
form of the gravitationally induced concentration gradient. If, as 
we suggest: 


Exponent 

Y 

6 

Of 

Relations: 








and dX/dZ parallels this behavior then: 

(X-X c )°oc(^ - z c ) (3.38) 

where z c would be the height at which the tricritical concentra¬ 
tion occurs. This implies that a careful determination of the 
functional relation between (ST — X c ) and (z — z c ) would not only 
give us values for (dX/5A) TJ> but would allow us to determine the 
critical exponent 6 as well. 

We will now consider the detailed prediction of several model 
calculations. Interest will be focused on how well the model 
reproduces the phase diagram of He 3 -He 4 mixtures, the extent 
to which behavior induced by changes in pressure is portrayed 
correctly, and the value provided by the model for the critical 
exponent 6. 

C. The Hard-Sphere Fermi-Bose Mixture 

The first model we will consider is that of a mixture of 
fermions and bosons with only hard-sphere interactions between 
the particles. This model was developed by van Leeuwen and 
Cohen 22 in 1968 and further elucidated by Kincaid and Cohen 23 in 
1975. Specifically they considered an anisotopic binary mixture 
of N b hard-spheres of zero spin with mass m B and diameter a 
mixed with N p hard-spheres of spin \ with mass m F and diameter 









a. This approximates the interatomic repulsive interactions by 
a hard- sphere potential and ignores any attractive interactions 
between helium atoms. As a result of these approximations it 
was later found necessary to apply a temperature-dependent ex¬ 
ternal pressure to the system so that its density would remain 
comparable to that of an isotopic helium mixture in similar cir¬ 
cumstances . 

The Hamiltonian considered was; 

- ( 3 - 39 ) 


with; 
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(3,40) 


and 



Where H u is simply the Hamiltonian for a mixture of an ideal 
Bose and an ideal Fermi gas, and H 1 contains the various inter¬ 
action energies between the bosons and the feimions. Since we 
are considering only hard sphere interactions, the various <p’S 

in H all have the property 







(3.42) 


^ = 0 if |r x -r a |>a 
tp = » if ]r x - r z \<a 


This Hamiltonian is solved using the pseudopotential method 
of Huang, Yang, and Lee 24 , This involves the substitution of 
pseudopotentials of the form 





(3.43) 


into H L . These forms were developed first by considering the 
two-body case and constructing a set of pseudopotentials that 
gave the correct energy levels. An W-body pseudopotential was 
then built up from the two-body pseudopotentials. The JV-body 
pseudopotential, which has the sphere diameter u a” as an ex¬ 
pansion parameter, is given only to first order in u a*\ 

The Hamiltonian is then solved by viewing H L as a perturba¬ 
tion on H ot and obtaining the energy levels, E ni correct to order 
“a". From these energy levels the partition function, and hence 
the free energy are calculated. Appropriate derivatives of the 
free energy, F, then give the quantities of interest: 





(3.44) 
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The phase diagram is then calculated by identifying the con¬ 
densation of the Bose component of the mixture with the onset of 
the superfluid transition in He 3 -He 4 mixtures and mapping out 
the two phase region by setting the chemical potentials of the 
respective components in the upper and lower phases equal. As 
mentioned earlier, in order to compensate for the lack of attrac¬ 
tive interactions, an external pressure has to be provided. A 
temperature-dependent pressure was chosen such that the 
density of the Bose component of the mixture was the same as 
that of He 4 liquid at OK. 

The results of these calculations appear in Fig. 5. The quali¬ 
tative agreement between the experimentally determined values 
(doited lines) and the theoretical calculations (solid lines) is 
remarkable, considering the relative simplicity of the assump¬ 
tions of this model. The predictions of the model for the effects 
of pressure variation are, however, not good at all. They are 
summarized in the table below and compared with experimental 
determ inations: 

Table 3.2 


As P increases: 

Hard Sphere Model 
T x increases 
T c /T x increases 


He 3 -He 4 mixture 

T x decreases 

T /T, first decreases, 
c then increases 


X, at r=0K decreases 


X t at 


T = OK first increases, 
then decreases. 


!N 
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Figure 5 


Phase diagram of the: Hard-Sphere Fermi-Bose Mixture model 
together with the SVP helium phase diagram 
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This lack of agreement is hardly surprising considering that all 
attractive interactions, which might be expected to affect the 
behavior of the sample as the average interatomic distance 
varied, were neglected. 

This model also supplies us with a value for the critical ex¬ 
ponent, 6. This fairly detailed calculation was carried out by 
Kincaid, 25 and gives two values for 6, one for X>X C (6 + ) and one 
for X<X C {6'): 

6+ = 1 {3.45) 

6 ‘ =2 

It can be shown that this variation in the value of 0 is the re¬ 
sult of certain assumptions necessary to the solution of the 
Hamiltonian. More careful expansions of the appropriate 
thermodynamic variables, as we shall see later, give 6 an un- 
equivocal value of 2 for all values of X near the TCP. 

D. Magnetic Models 

The first magnetic model we will consider is that of Elume, 
Emery and Griffiths 36 developed in 1971. They model He 3 -He 4 
mixtures by considering a lattice, each of whose sites is oc¬ 
cupied by a fictitious spin variable S t with possible values 0 and 
± l. a He 3 atom at a particular site corresponds to the =0 

state and a lie 4 atom to S,dbl- All sites are occupied by one atom 
of one type or the other. The spin associated with He 4 is double 
valued to allow for the possibility of superfluid ordering. In fact, 


1 
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the superfluid order parameter is defined as: 

1 * 

vZ ^ O = thermal average (3.46) 

JV 1 Pl 

With M -0 indicative of normal fluid and M ± 0 indicative of 
superfluid. The number of He 3 and He 4 atoms present is given 
by: 


d-S( 2 > 


I -i 




and the He 3 concentration by; 


(3.47) 


X* 


m 

N 


(3.48) 


If the system is translationally invariant, this gives: 
M= <S ( ) 


(3.49) 


and 

X-l ™(S, a ) ( 3 - 5 °) 

The Hamiltonian considered is: 

H=H S +H t ~ t3 - 51) 

where and /i 4 are the chemical potentials of He 3 and He 4 and 


sfoastf .{UGIS 








35 


N s J 2D S * S i {3.52) 

JTft 1 

tin = summation over nearest neighbor pairs 

Since this term in the Hamiltonian will lead to a second-order 
phase transition at which M spontaneously assumes a non-zero 
value, it is responsible for the "superfluid ordering” of the 
sample. In order to include the possibility of other sorts of in¬ 
teractions being responsible for the behavior of their model, the 
authors include an interaction term; 

m mi 

+ (3.53) 

nn 

supplying different interaction energies if neighboring sites are 
occupied by He 3 -He 3 , He 3 -He 4 , or He 4 -He 4 pairs. 

The total Hamiltonian is then rewritten; 

H=- jJ2 s r s i~ # 2D $1S/ + A 2D S * “ N ^ zK ™ + ^ ( 3 - 54 ) 

nn nn i 

where: 

f? =^33 + ^44 - 2K 34 
A = p 3 - ii 4 + 2z (K 33 - K 34 ) 

and; 

z - number of nearest neighbors , 
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By varying the relative strengths of J and K, the types of 
phase diagram generated by an entirely "superfluid driven" 
system {K = 0, t7 # 0) , or that of a system whose behavior is en¬ 
tirely determined by the relative strengths of interaction between 
its constituents {J = 0 f K4 0), or that of any intermediate case, 
can be explored. The authors solve this Hamiltonian in the 
molecular field approximation and arrive at the following results 
for X and Af; 


, _ 2 coshj35M 

1 " exp{ j3[ A -tK{ 1 - X) ]} + 2 eosh3 3M 

2 Sinh(3 M 

_ exp{/3[A -K(l -X)]} -f*2cosh;3$M 


(3.55) 

(3.56) 


3=Jz 

j3=(fer) _1 


In general, it is necessary to solve these self-consistent equa¬ 
tions numerically, assuming values for 3, X, and A, There may 

f 

be several pairs of answers (M, X) for any given set of param¬ 
eters, and those are chosen, which minimize the free energy. 
Employing these methods, Blume, Emery and Griffiths exploicd 
the phase diagram for several ratios of £f/X. Fortunately, the 
phase diagram that most resembles that of He 3 -He mixtures oc¬ 
curs whenXsO, thus uncoupling eqs. (3.55) and (3.56) and 
making their solution that much easier. A phase diagram with 
X=0 is displayed in Fig, 6. A region of normal fluid with A/=0, 
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Figure 6 


The Blume, Emery and Griffiths 26 model phase diagram when 
#=0. The reduced temperature axis is in units of T/3. 
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is separated from a region of superfluid (M# 0) by a lambda line, 

designated 7 s (x). In this model, the lambda line is straight, and 
given by: 

T s =3 (I-* - ) (3.57) 

For 7/,3 values below 1/3, two sets of solutions to eqs. (3.55) 
and (3,56) minimize the free energy and the mixture separates 
into two phases, one, with M = 0, designated by X^X 0 {T) and 
one, with M*0, designated by X-X^T). The coexistence curve 
is topped by a sharp angle at A, the TCP, and the lambda line 
intersects the coexistence curve exactly at A, The coordinates 
of A are T c /3 = 1/3, X c = 2/3. 

Obviously, the model displays very good qualitative agree¬ 
ment with experiment. Even its quantitative agreement is quite 
good as the trueX c is 0.669 ± .005 and T c /T x is 0.402± .005. The 
main differences between model and experiment are: the lack of 
any T = OK solubility for He 3 in He 1 because the Fermi statistics 
of He 3 were ignored, and the continuity of slope between the 
lambda line and the right-hand branch of the coexistence curve. 
This latter has been proved an artifact of the molecular field ap¬ 
proximation and a careful series expansion of the parameters of 
this model by Saul and Wortis 27 results in differing slopes for 
the lambda line and the right-hand side of the coexistence curve. 

This model is not well suited for any determination of its 
pressure sensitivity. As all lattice sites are filled, no variations 
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in density are permitted, and no terms involving pressure or 
volume appear in the thermodynamic potentials. We will there¬ 
fore later discuss a similar model due to Takagi that does permit 
lattice vacancies, and therefore a determination of its behavior 
at various densities. 

We must finally make a determination of the value for the 
critical exponent, 6, predicted by this model. This is done ex¬ 
plicitly from eqs. (3.55) and (3.56) in Appendix C with the result 
that 6* = 1, 5~ = 2, just as in the hard sphere model. These cal¬ 
culations also result in a functional form for A in terms of X 
near the TCP. We are thus able to calculate (aA/ax} TP for 
various values of the adjustable parameters $ and K. Therefore, 
even though the calculation of an explicit pressure dependence is 
impossible in this model, it is tempting to adjust the above pa¬ 
rameters so that T c and X c take on the values appropriate to 
different pressures. We might then hope that the concentration 
susceptibility calculated from these parameters will give some 
idea of the pressure dependence of the gravitational gradient. 

The result of such a procedure is that at a given value of (X-X c ), 
(3 x/9A) tp increases as X c increases. Remembering the experi¬ 
mental behavior of X c as a function of pressure, this indicates 
that the gravitationally induced concentration gradient should be 
a maximum at P ^9.9 atm and decrease monotonically as the 
pressure moves away from that value. We will later discuss 
this prediction in the light of our experimental determinations. 
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The other magnetic model we will consider is that of Takagi 28 . 
Outwardly similar to the Blume, Emery and Griffiths model, it 
is a generalization of the work of Matsubara and Matsuda 29 , who 
dealt with a lattice model for pure He 4 . In order to incorporate 
the dynamics of the He 3 without increasing the degrees of free¬ 
dom of his model, Takagi treats the He 3 as a classical gas 
"background”, solves its dynamics separately, and includes its 
kinetic energy term as part of its chemical potential. This lack 
of parallelism in the treatment of He 3 and He 4 will cause us 
severe difficulties later, when we try to extract a value of A 
from this model. 

In this model, again, we have a lattice, each of whose sites is 
occupied by a fictitious spin. The z components of these spins 
have the range of values: -1,0,1. The assignment of physical 
significance to these spin values is, however, very different than 
in the Blume, Emery and Griffiths model: +1 indicates a site 
occupied by He 4 ; 0, a vacant site; and -1 a site occupied by He 3 . 
This is not an Ising model in that the possibility of spin ordering 
in the X-Y plane is not excluded. Indeed, it is just such an 
ordering that produces the superfluid transition. This assign¬ 
ment of spins leads, not surprisingly, to a very different 
Hamiltonian: 

H=H S +H I +H C ( 3 - 58 > 


with: 



(3.59) 



( 3 * 60) 


jm 



(3.61) 


Where y 0 is the number of nearest neighbors, K, I, and L are 
determined by the values of the various intraparticle potentials, 
and H l and H 2 are functions of the chemical potentials including 
a kinetic energy term for the He 3 . H s is a superfluid ordering 


term resulting from the dynamics of He*. H T is an interaction 


term whose coefficients are combinations of the potentials V u 
between He 1 and He J . H c is a term that includes the effects of 
the chemical potentials of He 3 and He 4 . This Hamiltonian is first 
presented in second quantization formalism and then cast in the 
above form by demonstrating that the appropriate combinations 
of spin operators are, in fact, the creation, destruction, and 
number operators required. Unfortunately, before going on to 
solve this Hamiltonian, Takagi sets H r = 0 and thus consideis 
only superfluid driven effects on the phase diagiam . 

The results of these calculations are shown in Fig. 7, for 
various densities. Again, the qualitative features of the phase 
diagram are duplicated quite well. We again observe the con¬ 
tinuity of slope between the lambda line and the right-hand 
branch of the coexistence curve, and the lack of solubility for 





Figure 7 


The phase diagrams of the Takagi 28 model for different lattice 
densities, p. Curves 1, 2 and 3 refer to p =0.65, 0.75 and 0.85, 
respectively. The reduced temperature axis is in units of T/ 
7\(X>0,p=i) 
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He 5 in He 4 at T= OK that occur for the same reasons as advanced 
for the Blume Emery and Griffiths model. Takagi also correctly 
predicts the pressure dependence of the lambda line-but this in 
itself is not too surprising, as the work of Matsubara and Mat- 
suda had already displayed the correct pressure dependence for 
T x of pure He 4 . Thus, this model only had to provide for a 
depression in T x with increasing X to be correct. Of more im¬ 
portance is this modePs prediction of the motion of the TCP with 
increasing density. Here , the prediction of a T c at first in¬ 
creasing and then decreasing with increased pressure while X c 
moves monotonically to lower values is not at all borne out by 
experiment. Quantatively also, this model is less effective than 
that of Blume, Emery and Griffiths. The coexistence curves 
are shaped less like the experimental ones, andJ£ c only ap¬ 
proaches its correct range of values as p — 0.95, where p is the 
number density of occupied sites. 

It was our hope to extract a value of 6 from this model and 
perhaps even a measure of how the strength of the gravitational 
gradient might vary with pressure. Proceeding analogously with 
our treatment of the Blume, Emery and Griffiths results, we 
were plagued by the lack of symmetry in the treatment of the 
chemical potentials. We did manage to extract a value of 6-2 
for approaching the superfluid side of that TCP foi which 
= \ j along a line on which p. A was kept equal to § J by varying 


both p andX. 
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Siemann and Huber 30 have investigated the critical behavior of 
the Takagi model near the lambda line and the TCP, Employing 
a method of double-time Green’s functions, they extract 
"classical” values for all of the critical exponents, in particular, 
6 = 2 in the "paramagnetic” (i.e. normal fluid) phase at the TCP 
along a line normal to the phase separation curve in T- A space. 
No determination is made of exponents in the superfluid phase, 
and the pressure dependence of (8 x/&A) tp is not investigated. 

In order to remove the apparent discontinuity in 6 predicted 
by the theories examined here, we refer to a paper by Riedel 
and Wegner 31 who employ methods developed by Wilson 33 to 
closely examine the region near the X line and critical point of 
the Blume, Emery and Griffiths model. They believe that the 
major cause of such discrepancies is the difficulty in achieving 
discontinuous behavior in the derivatives of the thermodynamic 
potential, when the potential itself is an analytic function. They 
claim that the usual methods employed, going to the “thermo¬ 
dynamic limit" of infinite volume, while correct, makes it very 
difficult to properly extract the singular behavior of the system 
and hence the critical exponents. To overcome this, they employ 
a generalization of the "block” scaling picture developed by 
Kadanoff 33 , extracting differential equations for the coupling 
constants and magnetic field appropriate to the block which 
lead, in a natural way, to the lack of analyticity of the partition 
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The process, as applied by Riedel and Wegner, involves a 
complex division and sectional integration over the phase space 
of the model and results in a constant value of 2 for 5 on both 
sides of the TCP. 

E. The Form of the Gravitationally Induced Concentration 

Gradient 

Assuming that 5 does indeed have the value 2, as both theory 
and experiment indicate, we can now show how the concentration 
gradient ought to behave in a real He 3 -He 4 mixture as a function 
of temperature. Fig. 8. Assume that the average concentration 

5g 

of the mixture is X c . Then for T» T c no appreciable gradient 
will be present. As T—T c , the gradient will begin to form, with 

a point of inflection atX c . When T - T c> parabolas with horizontal 

* 

tangents at X~X C will form in both the upper and lower phases, 
with their common origin marking the incipient phase boundary. 

Finally, at T< T Cf the phase boundary is formed, and gravita¬ 
tional gradients exist in both the upper and lower phases. Note 
that now neither parabola displays a horizontal tangent as (9A/ 

9x) tp =0 only at the TCP. The gradient in the He 3 rich phase is 
shown as smaller than that in the He 4 rich phase since there is 
experimental evidence 12 ' 13 * 18,19 that it is less by about an order of 
magnitude. Our measurements have confirmed this. It should be 
mentioned that the behavior of a real He 3 -He 4 mixture is more 
complicated than that portrayed in Fig. 8, because the average 






Figure 8 


A qualitative representation of the gravitationally induced con¬ 
centration gradient in He 3 -He" 5 mixtures at several temperatures 
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concentration of such a mixture is constant. Thus the height of 
the interface and the "incipient interface” is not constant. It is, 
in fact, a very strong function of the temperature, the magnitude 
of the concentration gradient, and the average concentration. 

In Fig. 9 we display the motion of the interface for T < T c 
using the concentration gradient parameters we have measured 
at 2.5 atm. Curves are show for various values of X-X c , It 
is the great sensitivity of the location of the interface to changes 
X-X c near the TCP that allows us to position it advantageously 
in our apparatus. The curves were generated by calculating the 
average value of X in the sample chamber given the concentration 
profile determined by the coexistence curve and the concentra¬ 
tion gradient. The height of the interface was then adjusted until 
X took on a specified value for a given T. Figure 9 is simply a 
plot of these height vs. temperature curves. 
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Figure 9 


Interfacial motion in He 3 -He 4 mixtures as a function of T- T c 
for various values of X -X c . 
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IV. Apparatus 
A. General 

The design of this apparatus was an attempt to meet the 
various, often conflicting, requirements of this experiment. The 
requirements were: 

1. To obtain a sometimes superfluid sample at pres¬ 
sures up to 30 atmospheres, 

2. To allow a determination of the pressure of the 
sample to within 0,1 atmospheres. 

3. To permit variations in the pressure of the sam¬ 
ple without changing its concentration. 

4. To maintain the sample at temperatures between 
500 mK and 900 mK for periods up to 24 hours. 

5. To measure these temperatures to within a few pK. 

6. To stabilize temperature, when required, to with¬ 
in a few tiK. & 

7. To minimize thermal gradients within the sample 
chamber. 

8. To determine the concentration of the sample at 
a particular height in the sample chamber, rela¬ 
tive to its average concentration, to within a few 
hundredths of a percent. 

9. To be able to vary and measure the vertical dis¬ 
tance between the concentration sensor and the 
phase-boundary of the sample. 

These requirements were met with varying degrees oi success 
by a sample chamber of quite complex design mounted in a stan¬ 
dard He 3 refrigerator cryostat. We decided to measure the 
concentration by the capacitative technique developed by Heybey 
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in 1962 34 employing a Tunnel Diode Oscillator (TDO). In this 
way we hoped to achieve the required sensitivity with a rela¬ 
tively simple design. To meet requirement #8, a single 
parallel plate capacitor was used as the concent ration-sensing 
element. The capacitor plates had to be horizontal and as closely 
spaced as possible. 

Requirement # 9 necessitated the most elaborate design. The 
capacitor measuring sample concentration had to remain sta¬ 
tionary for greatest stability; therefore we contrived to move 
the phase boundary instead, by installing a movable plunger in¬ 
side the sample chamber. By moving this plunger up and down, 
the location of the phase boundary could be vertically shifted 

■ 

relative to the location of the capacitor, as shown in Fig. 10. 

We monitored temperatures by using semiconductor resistance 
thermometry, measuring these resistances with AC bridge net¬ 
works. Temperatures were controlled by allowing a suitable 
bridge output to drive a current source that heated resistance 

wires surrounding the sample chamber. 

The design of each subsystem of the apparatus is detailed in 

the following sections. 

B. Cryogenics 

Low temperatures were maintained by two nested dewar s of 
silvered glass manufactured by Kontes/Martin. The outer dewar 
ol 5.5" IJ). was filled with liquid nitrogen while the inner dewar 










Figure 10 


An illustration of how plunger motion can affect the relative 
positions of the interface and the concentration sensor. 
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of 3" LD. contained a liquid helium bath. This inner dewar was 
an integral part of a vacuum system which allowed the bath to be 
pumped to an ultimate temperature of *'1.1K. This temperature 
could be maintained for at least 24 hours. Other components of 
this vacuum system included a combination manometer-over¬ 
pressure relief value , a 150 CFM Kinney KT-150 vacuum pump 
connected to the dewar by 4" pyrex lines, and a 5 CFM Direc- 
Torr 8805 vacuum pump that was used to maintain the bath at 2 
K and higher when data was not being collected, A fail safe sys¬ 
tem was devised to protect the Kinney pump from interruptions 
in the supply of electricity and cooling water and the dewars 
from failures of the two vacuum pumps. This allowed data runs 
to continue without the constant attention of an observer. 

C. Cryostat 

Suspended in the helium bath by two T- 304 stainless steel 
pumping lines was the cryostat. Fig. 11. Its main components 
were an exchange-gas can 2.5" in diameter and 10" long, de¬ 
mountable via a Wood’s metal joint, and an oxygen-free high 
conductivity (OFIIC) copper platform with an integrally machined 
He 3 reservior of about 4 cc capacity. Pumping on the He 3 in this 
reservoir permitted the maintenance of the copper platform and 
any attached equipment in the range from 1 K to 500 mK for 
periods up to 24 hours. The sample chamber described below 
was attached to this platform with varying degrees of thermal 










Figure 11 


The Cryostat 
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contact during the course of the experiment. Thermal grounds 
were provided both on the exchange gas can (at bath temperature) 
and on the copper platform {at He 3 temperature). 

The top of the exchange gas can was pierced by a variety of 
vacuum-tight feed-throughs to allow access to the experiment. 
Most of the elctrical connections were made via advance wires 
that ran down inside the exchange gas pumping line. The only 
exception to this was a coaxial cable feed-through made-as were 
all other electrical feed throughs-by passing forravar insulated 
copper magnet wire through small diameter (£" or less) copper 
tubing and then vacuum-impregnating the assembly with filled 
epoxy (Emerson and Cummings 2860 FT). The body of a low- 
temperature valve also passed through the exchange gas can as 
did a parallel sample lines of .004" I.D. stainless steel capillary. 
A rotating shaft ran down the central axis of the cryostat and 
passed through the exchange gas can and the copper platform. 

The vacuum seal for this shaft was at the top of the cryostat at 
room temperature. 

The two pumping lines mentioned above were connected to two 
separate vacuum systems. The exchange gas line was connected 
through a liquid nitrogen cold trap to a 4" Veeco 41W diffusion 
pump backed by a Duo-Seal Model 1400 mechanical pump. 
Monitored by both a thermocouple and an ion vacuum gauge 
this system commonly achieved vacuums of 5 X 10- Terr in the 
exchange gas can. The lie 5 system was quite similar except that 
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the exhaust of the mechanical pump-an especially sealed duo- 
seal model 1402-ms collected in two 1.5/ metal storage 
vessels for re-condensation into the He 3 reservoir. This sys¬ 
tem also incorporated a Hoke 2331 needle value between the 
diffusion pump (which was actually never used) and the He 3 
reservoir. By means of this valve the rate of evaporation of 
He 3 could be controlled to give cooling rates in the neighborhood 
of I mK/Hr in the required temperature range. 

D. Sample Chamber 

The sample chamber, Fig. 12, ms mounted on the copper 
platform by a threaded brass ring Thermal contact between 
platform and chamber was provided by several adjustable 
copper straps. All electrical leads were first thermally 
grounded at the bath temperature and then to the He 3 reservoir 
before entering the sample chamber. The central shaft, of 1/16 
stainless steel, was extended by 4" of Micarta shafting before 
being attached to the sample chamber. These steps were taken 
to reduce the heat input to the sample chamber by these routes 
to a minimum, thus reducing the possibility of thermal gradients. 

The body of the sample chamber was machined from a 2.25 
diameter cylinder of OFHC copper. The top section of the sam¬ 
ple chamber contained the magnet drive, whereby the rotary 
motion of the central shaft was transmitted inside the sample 
chamber. The relatively large diameter of the chamber was a 
constraint imposed by the diameter of reasonably strong magents. 
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Figure 12 


A cross sectional view of the sample chamber 
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The magnets employed - alnico disks of 1.63" diameter and .26" 
thickness- were capable of transmitting a torque of several 
inch-ounces across the 1/16" stainless steel diaphragm sepa¬ 
rating them. Both the upper and lower magnets were supported 
by ball bearings. It was necessary to increase the clearances 
of these bearings to insure smooth running at liquid helium 
temperatures. This was done by running them at room tempera¬ 
ture while filled with a mixture of fine grinding compound 
(600 grit) and oil for a minute or two. The magnet drive rotated 
a 4-40 stainless steel threaded rod that ran the length of the sam¬ 
ple chamber and was supported by an adjustable needle bearing at 
the other end. 

Threaded on this rod was a plunger assembly that was keyed 
into a slot in the sample chamber wall. The cross sectional area 
of this plunger was one half the free cross sectional area of the 
bore of the sample chamber. Thus the rotation of the central 
shaft caused a rotation of the threaded rod which, in turn, moved 
the plunger up and down along the bore of the sample chamber. 

The ratio of the areas was chosen to make the motion of the inter¬ 
face equal and opposite to the motion of the plunger. The length 
of the plunger was made f the length of the bore of the sample 
chamber to maximize the distance through which the interface 
could be moved. The maximum travel of the plunger was ■% " 
i.e. 20 turns of the drive screw. 
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Hie capacitor sub-assembly, used for making dielectric mea¬ 
surements on the sample, was designed to contain sample over 
as small a vertical height as possible, while attempting to keep 
stray capacitance to a minimum. It was manufactured by 
fastening the upper capacitor plate, with its attached electrical 
connection, into Us proper position on the base plate (maintained 
by a 0.002" plastic shim) ’with sealing wax. The "epoxy fill” 
(Emmerson and Cummings 2850 FT) was then added while the 
entire assembly was enclosed in an appropriately shaped mold 
coated with a release compound. Warming the capacitor assembly 
after the epoxy had set allowed removal of the shim, and acetone 
removed the last traces of sealing wax. This process resulted 
in a parallel plate capacitor of ~ 0.004" spacing and with a fill 
factor (the ratio of active to stray capacitance) of better than j . 
The capacitor sub- assembly was then mounted in a recess that 
communicated with the bore of the sample chamber. This design 
was chosen to minimize pressure effects on the capacitor. Once 
connected in parallel with a toroidal coil of tinned copper wire 
wound on a Micarta form, the capacitor became part of a tank 
circuit with a resonant frequency of ~11.5 MHz. 

A removable bottom plate on the sample chamber gave access 
to the bore and plunger area, and contained a germanium resistor 
(in a He 3 -filled capsule) and its associated electrical feed- 
throughs as a temperature sensor. All the various sections of 
the sample chamber were assembled with 0.005" lead O- rings 
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and brass 4-40 fiilister~head screws to create an hermetically 
sealed unit with an interior sample volume of about 3 cc. 

E. Electrical 

The tank circuit mentioned above was part of a tunnel diode 
oscillator (TDO), Fig- 13. The circuit elements of the oscillator 
were divided into two sections, the biasing resistors being 
mounted on the exchange gas can heat sink, as they were the 
main dissipators of heat (100 4 W), and the diode and associated 
passive elements being mounted directly on the sample chamber. 
In order to keep the heat input to the sample chamber as small 
as possible, the diode was attached to the tank circuit by ~-l cm 
of #39 advance wire, while it was connected to its voltage supply 
by four strands of # 18 copper wire thermally anchored to the 
He 3 pot. It was thus our hope that the majority of the ~2 pw 
dissipated in the tunnel diode would not enter the sample chamber 
at all. 

The TDO was powered by a 1.35 V "D" size mercury battery. 
In order to minimize temperature effects upon this cell it was 
mounted on a 500 g aluminum block and surrounded by one inch 
styrofoam walls to create a long thermal time constant. This 
was necessary because the TDO frequency varies slightly ac¬ 
cording to the biasing voltage on the diode. The output of this 
cell could be fed either directly to the biasing resistors of the 
TDO or through a 25 KO potentiometer, enabling us to adjust 










Figure 13 


Schematic diagram of the Tunnel Diode Oscillator 
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the biasing voltage. This biasing adjustment was necessary only 
when we ran the TDO at room temperature for adjustment pur¬ 
poses. 

The signal from the tunnel diode was carried to the top of the 
cryostat by the coaxial line. This was made from a piece of j *' 
O.D. 0.01" wall, 304 stainless steel tubing containing a teflon 
tubing insulator and a center conductor of alternating copper and 
advance wires. The top and bottom 6" of the center conductor 
were of advance. At the cryostat head, the signal was fed into a 
tuned amplifier, Fig. 14 and thence into a wide band amplifier - 
an Hewlett Packard, (H.P.), 8447 A. This resulted in a 0.3 V 
rms signal that was counted by an H.P. 5245 L counter. The 
counter output was then fed into a digital-to analog converter 
(H.P. 580A) and displayed on a 10" strip chart recorder 
(H.P. 7100 BM). 

The TDO ran at -11.5 MHz at 1.1 K with a stability of ±0.5 Hz 
over several hours with the sample chamber empty. The depen¬ 
dence of frequency upon temperature was of the same order - 
0.12 Hz/mK. Since a typical run would cover a range of about 
10 raK, this would lead to frequency corrections on the order of 
1 Hz. Long term stability was affected by sudden frequency 
jumps of ”30 Hz which occurred only during the thermal cycling 
associated with refilling the He 4 bath. Thus, the intrinsic 
stability of the TDO during a run was about ± 1 Hz, not including 
the effects of actual concentration instabilities in the sample. 
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Figure 14 


Schematic diagram of the tuned amplifier 
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Between runs, frequency changes of "-SO Hz were possible. 

F. Temperature Measurement 

Three temperature measuring elements were attached to the 
cryostat and sample chamber. A Speer 100 $l, \ watt carbon 
composition resistor, grade 1002, was attached to the exchange 
gas can. A similar resistor was mounted on the He 3 reservoir, 
and a Germanium resistor (Cryocal CR500-He 3 ) was mounted 
inside the sample chamber. 

Two A.C.bridge networks were used to measure these re¬ 
sistances. One, of relatively low sensitivity was switchable 
between the two Speer resistors, and one of high sensitvity and 
stability, Fig. 15, measured the resistance of the Germanium 
thermometer. The effect of changes in room temperature upon 
the calibration of this bridge were kept to a minimum by em¬ 
ploying a fixed, known resistor of low temperature coefficient 
(a Vishay type S102 with a coefficient of 0±1 ppm/°C between 
0 - 60 °c) and adjusting the balance of the bridge with a ratio 
transformer (Gertsch model 1004) with a negligible temperature 
coefficient. Any possible change of bridge calibration due to a 
change in the resistance of the wires Leading to the Germanium 
resistor was compensated for by a set of identical leads included 
in the circuit with the known resistor. An unanticipated benefit 
of this sort of circuit is that readings given by the ratio trans¬ 
former were approximately linear in temperature over a very 









Figure 15 


Schematic diagram of the A.C. resistance bridge circuit 
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wide range (600 - 900 mK), even though the resistance of the 
Germanium thermometer changes exponentially with tempera¬ 
ture. 

To eliminate ground loop noise and other types of interference, 
the entire bridge circuit was kept floating by the use of input and 
output transformers. All leads were twisted pairs run inside 
grounded shields. The bridge circuits themselves were built in¬ 
side a grounded steel chasis. The signal from the bridge circuit 
was fed back to the PAR 124 lock-in amplifier through a balanced 
input to the model 119 preamp, in the "direct” mode. Impedances 
were matched by the output isolation transformer to give the best 
noise figure. 

The Speer resistors were calibrated only for the temperatures 
of boiling liquid nitrogen (77 K) and liquid He 4 (4.2 K), for the 
superfluid transition of He 4 (2.17 K) and for the minimum bath 
temperature (-1.1 K). They were mainly employed for coarse 
determinations of the temperature and level of the He 4 bath during 
transfer of liquid He 4 . The Germanium resistor was calibrated 
against the vapor pressure of pure He 3 as given in the T 6 ' tem¬ 
perature scale of Roberts et al-? s and thereafter served as the 
primary temperature scale during the course of the experiment. 
This calibration was accomplished by measuring the pressure 
over the He 3 reservoir with an M-K.S. Barytron pressure meter, 
type 1200, equipped with a model 145 AH head capable of mea¬ 
suring 10 mm Hg full scale with an absolute accuracy of a part 
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in IQ 4 . The Baryfcron itself was calibrated both against an oil 
manometer and a dead-weight pressure gauge manufactured by 
the Ruska Instrument Corporation, whose accuracy is traceable 
to the National Bureau of Standards. 

The Barytron, besides giving a most accurate and convenient 
method of measuring the He 3 vapor pressure, allowed us to con¬ 
stantly observe and record this pressure on the strip chart 
record of the experiment, thus guarding against changes in the 
calibration of the Germanium resistor and associated circuitry. 

Temperature calibrations were carried out in both slowly 
warming and cooling modes at rates not exceeding 4 mK/Hr. 
These results were averaged to compensate for any temperature 
differential between the sample chamber and the He 3 reservoir. 
The sensitivity of the bridge allowed relative temperature mea¬ 
surements of 4 pK and the calibration proceedure resulted in an 
accuracy of ± .3 mK relative to the T 62 scale. 

The final element of the temperature measuring system was 
the "zero finder”. This was simply an operational amplifier 
driven by the function output of the PAR 124. Because of the 
extremely large amplification factor of the operational amplifier, 
its output would essentially always be driven to one of its ex¬ 
tremes, ±10 V, except when the input signal passed through 
zero, at which point the amplifier output would jump from one 
extreme to the other. This output voltage was channeled by 
diodes to operate a relay when the amplifier output went negative. 


H J*V #** . 







The i elay had both normally-open and normally-closed contacts 
which were used to energize the event marker on the strip chart 
recorder, and a buzzer warning the experimentalist that the 
desired temperature had been reached. 

G. Temperature Control 

Temperature control was accomplished by allowing the output 
of the PAR 124 to control a current source, whose output could 
be divided among any of four heaters inside the cryostat. All 
healers were wound of #39 advance wire, one on the He 3 reser¬ 
voir and three at various levels on the sample chamber. By thus 
distributing the heat input to the sample, we hoped to minimize 
the effect of any thermal gradients thus induced. This system, 
when driven by the bridge at a PAR 124 sensitivity of 2 p.V with 
a 3 second time constant could maintain temperatures to within 
20 pK, A manual heating feature allowed coarse temperature 
adjustments to be made by hand. 

H. Pressure Measurement and Control 

Pressures were both maintained and measured by a high 
pressure gas handling system, Fig. 16. Built entirely of stain¬ 
less steel, the system used Nupro SS-4H bellows valves con 
connected to }" tubing of 0.030" wall by Swagelok joints. The 
reservoirs were 1 liter high pressure sample cyclinders, rated 
at 1500 pSlG, with i " NPT connectors at each end. 
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Schematic diagram of the high pressure gas handling system 
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The system used nitrogen gas to force a mercury “piston” 
to compress the sample to the desired pressure. Heise Bourdon 
tube gauges model C 58150 with a range of 0- 30 atmospheres 
and a sensitivity of 0.01 atm. were used to measure the pressure 
on either side of the mercury. The gauge that measured sample 
pressure was calibrated against the above mentioned dead weight 
gauge and no deviations in accuracy were found to within its 
sensitivity. The differential pressure between the two gauges 
was used to determine the mercury level within the system. 

I. Sample Lines 

AIL room temperature sample lines were of copper capillary 
tubing of 0.026" I.D. for low dead volume. Where necessary, 
they were connected to stainless steel Nupro SS-4H bellows 
valves by adaptors that consisted of j " stainless steel tubing 
to be attached to the valve's Swagelok fittings hard soldered to 
i" brass tubing whose I.D. was drilled to accept the copper 
capillary for soft soldering. The sample was conducted down the 
cryostat via two routes: a 0.004" I.D. stainless steel capillary 
valved at room temperature, and a low temperature valve, Fig. 

17, whose valve stem was the sample channel. Once the low 
temperature valve was closed, the sample chamber and associated 
4 mil capillaries were isolated from the rest of the system, thus 
maintaining the concentration integrity of the sample. 
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V. Procedure and Results 


A. A typical run 

Starting with both the He" 1 and liquid nitrogen baths filled and 
the exchange gas can evacuated to ~ 5x 10* 0 Torr, gaseous He 3 
was admitted to the He 3 reservior via the pumping line. The He 4 
bath was then pumped, reducing its ultimate temperature to 
~U K and condensing the He 3 gas into its reservior. At this 
time the He 3 was acting as a thermal link between the sample 
chamber and the He 4 bath. In general, a period of about thirty 
minutes was necessary for thermal equilibrium to be established 
at the ultimate bath temperature. 

When the temperature of the cryostat stabilized at “1.1 K, a 
He 3 -He 4 mixture would be condensed into the sample chamber 
by way of the low temperature value. Pressures varying between 
1.5 and 6 atmospheres were used to force the sample through the 
somewhat constricted geometry of its filling channel to reduce 
the possibility o£ generating changes in concentration by fractional 
distillation. Typically, the sample chamber could be completely 
filled by this method in 15 minutes. Tho sample would then be 
pressurized to approximately the value desired and the low tem¬ 
perature valve closed. The top of the valve stem would then be 
isolated from the rest of the high-pressure system so that any 
subsequent fractional distillation of the sample within the valve 
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stem could not affect the concentration of the remaining room 
temperature mixture. 

The filling of the sample chamber, initially under vacuum, 
was graphically portrayed by the frequency readout on the strip 
chart recorder. As soon as the sample began to enter, the fre¬ 
quency would drop by about 3 KHz, indicating that the dielectric 
constant between the capacitor plates had been increased by the 
presence of gas and, perhaps, superfluid film. The frequency 
would then remain constant to within a few hundred Hz until the 
liquid level reached the height of the capacitor. The filling of the 
capacitor with liquid would cause a precipitous drop in frequency 
(150 KHz in 2 minutes). The frequency would then stabilize again. 
The completion of the filling operation would be marked by 
another sudden downward excursion in frequency as the sample 
became pressurized. This pressurization would be continued 
until an appropriate pressure range was reached. When the fre¬ 
quency finally stabilized again, it was concluded that no pressure 
gradients remained in the system and that the pressure read at 
the high pressure gas handling system accurately reflected the 
pressure inside the sample chamber. The closing of fhe low 
temperature valve would induce small frequency changes on the 
order of 20 Hz indicative of changes cither in pressure (-0.03 Atm) 
or, less likely, concentration induced by this action. 

The sample chamber would then be cooled by pumping on the 
He’ reservoir until a temperature of 870 mK reached. This 











temperature would be maintained by the temperature controller. 
The pressure was then trimmed by opening the valve on the 4 
mil capillary until a stable frequency again indicated pressure 
equilibrium. In general, the pressure of the sample would in¬ 
crease by about 0.3 atm. during this cool-down. The 870 mK 
temperature was chosen as appropriate for trimming the pressure 
as the sample could have no superfluid component at this tempera¬ 
ture over the range of pressures and concentrations considered. 

It was feared that initiating a flow of sample while some of it was 
superfluid would result in a "heat flush” effect that might drasti¬ 
cally alter the concentration. After again isolating the sample 
chamber from the high pressure gas handling system, it would 
be cooled to the range of temperature in which we anticipated 
that phase separation would take place. 

Our first aim was to determine the concentration of the sam¬ 
ple relative to the tricritical concentration at the 

pressure being studied. This was done by cooling the sample 
through the phase separation temperature twice, once with the 
plunger in its bottom position and once with the plunger at the 
top. if X «X C the onset of phase separation would be marked by 
a drop in oscillator frequency during both cooling runs. The 
drop in frequency was due to the increase in He- concentration of 
the lower phase. That both cooling runs showed similar drops 
in frequency was indicative of a phaso boundaiy coming into 
existence above the scanning range of the sample chamber. 










Analogously if X »X C the onset of phase separation would be 
marked by an increase in frequency during both cooling runs. 

Once the concentration of the sample relative to X* was 
determined, we would adjust the value of X-X c by changing the 
pressure of the sample by way of the 4 mil line. This would 
change X Q while, hopefully, leaving X constant. Our goal in this 
procedure was to finally arrive at a value of X-X c for which the 
two cooling runs described above would show opposite frequency 
effects, indicating that the phase boundary had formed within the 
scanning range of the sample chamber. Several cooling runs 
would then be done with the plunger in various positions in the 
sample chamber. Usually some small additional corrections to 
the pressure were necessary to cause the point of formation and 
subsequent motion of the interface with decreasing temperature 
to be appropriate for the taking of data. We would strive for a 
value of X-X„ that would form the interface relatively high in the 
sample chamber at T = T C and keep its motion to a minimum for 
temperatures slightly less than T c . In the case depicted in Fig* 

9, for example, we would attempt to achieve a relative concen¬ 
tration of -0.001. 

We had originally hoped that when the correct pressure had 
been achieved, it would be possible to Lake data by simply holding 
temperature constant and varying the plunger position, thus di 
rectly measuring the variation in concentration with the distance 
from the interface. Unfortunately, the heating and vibration 








associated with the rotation of the magnet drive were sufficient 
to induce permanent shifts of up to 30 Hz in the oscillator out¬ 
put. Since the effect for which we are searching has a breadth 
on the order of 100 Hz, it was clearly impossible to take data in 
this fashion. We therefore elected to take data by varying tem¬ 
perature while keeping the plunger position constant. We would 
begin at some temperature sufficiently high to guarantee that the 
gravitational effects were negligible and cool slowly at rates less 
than 4 mK/Hr. until a point about 10 mK below the phase separa¬ 
tion temperature was reached. The sample would then be 
warmed, the plunger moved to a new position, and the process 
repeated. This procedure allowed us to compensate for the fre¬ 
quency changes caused by the use of the magnet drive by adjusting 
our data so that the concentration measured at the highest tem¬ 
perature during the run, assumed equal to the average 
concentration of the sample, was always the same. To further 
guarantee that these adjustments were accurate, the concentra- 
tions measured by the capacitor below the temperature range in 
which gravitational elfects were Important also had to agree ior 
the data to be considered valid. Then data taken at the same 
temperature from runs at dUterent plunger positions could he 
correlated to explore the relationship between the desired 

variables. 











B. Methods of Data Analysis 

We have previously described the methods employed in 
determining the pressure and temperature of the sample. Rela¬ 
tive concentrations were measured by determining the concentration 
sensitivity of our TDO circuit. As is shown in Appendix D there 
is a linear relationship between changes in frequency and changes 
in concentration, eq. D-12. The constant of proportionality depends 
upon what pressure is being considered. The relative motion of 
the interface caused by movement of the plunger is a bit more 
complex. We assumed, as a first approximation, that the motions 
of the plunger and the interface were equal and opposite. This 
would be rigorously true only if the interface separated two phases 
of uniform concentration. From this assumption it was possible to 
determine an initial concentration profile of the sample at various 
temperatures. This concentration profile was used to calculate 
interfacial motion curves of the type shown in Fig, 9 for the dif¬ 
ferent internal geometries caused by different plunger positions. 
From a set of these curves at different plunger positions, the 
actual motion of the interface caused by changing the plunger posi¬ 
tion could be determined. The actual motion of the interface was 
discovered to be approximately linear in plunger position, but 
deviations in linearity by as much as 3% had a significant effect on 

the determination of critical exponents. 

Fig. 18 shows how typical data at different plunger positions 
appeared when plotted on T-X coordinates. These plots were 









Figure 18 


Typical data from cooling runs at different plunger positions 
plotted on T-X coordinates. The left hand branch of the co¬ 
existence curve is extended towards the TCP (dotted line) for 
the determination of T c and the relative value of X c , Only a few 
data points are shown, for clarity. 
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also used to determine T c and X c by extending the left and right 
hand branches of the coexistence curve towards the TCP. Fig, 
I9 f shows the same data plotted in z - X space where z is the 
vertical distance between the interface and some arbitrary 
reference height. Plots of this type were used to determine 

C J 

the height of the tricritical interface. The uncertainties in these 
graphical determinations were typically ±0.0005 in concentration 
and ±0.3 mm. in height. 

The data was fit to the functional form: 


A(X-X c Y = (z-z c ) 


(5.1) 


with A and 6 being determined by least squares analysis. X c and 
z c were allowed to vary within their -uncertainties and values 
chosen for them that resulted in a minimum for the sum of the 
squared deviations of the data. 8 was identified as the critical 
exponent of the system and A as a measure of the magnitude of 
the gravitational gradient . 


C. Experimental Results 


Data were taken at pressures of 2.5±0.1, 9.9±0.1, and 16.7±0.1 
atmospheres. While the experiment was being run at 2.5 atm. the 


sample chamber was in good thermal contact with the He 
reservoir. This caused temperature fluctuations in the sample 
that resulted In concentration fluctuations of ±0.001. The source 
of these temperature fluctuations seemed to be the He 3 reservoir 









Figure 19 


The same data as in Fig 18, plotted on X- z coordinates at 
different values o£ T. The height at which the interface disappear; 
determines z c ~ Interfaces are shown as dotted lines. Only a if’' 
data points are shown, for clarity. 
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itself. We assume that some defect in the metering valve or, 
perhaps, a superheating of the liquid He 3 in the reservoir and its 
subsequent relaxation by way of violent boiling resulted in an un¬ 
predictable cooling rate, In order to reduce these concentration 
fluctuations, the temperature was stabilized by means of the 
temperature control circuitry at each point desired. Unfortunate¬ 
ly, the activation of the heaters necessary for temperature control 
resulted in a consistent downward excursion of the concentration 
measured by the sensing capacitor of -0,001, The system did, 
however, become much more stable, showing concentration 
fluctuations of only ±0,002 while temperature was stabilized. We 
therefore attempted to usg data gathered in this manner, compen¬ 
sating for the deviation in concentration resulting from the 
activation of the heaters with an additive constant equal to the 
average value of these deviations. Analysis of these results led 
to a value of 6 *2.85±0.3 and to internal inconsistencies. It was 
found that values of X c and z c predicted by the best fit of the data 
to eq, 5.1 were well outside the perm is sable range of these 
parameters predicted by the graphical analysis of the data as 
described in section 5B. We therefore concluded that the use of 
heaters resulted in a heat flush effect, changing the concentration 
being measured by the TDO circuit in some complex and un¬ 
predictable way. The results were then reanalyzed, taking con¬ 
centration values at points where temperature was not being 
stabilized. The strip chart record was smoothed as best 






as possible by eye in order to compensate for the concentration 
fluctuations. Fortunately, the extrapolated value of was not 
one of the points where temperature stabilization was employed. 

As we shall see, this analysis gave results of greater internal 
consistency and in better agreement with previous work. 

The cryostat was modified by reducing the thermal contact 
between the sample chamber and the He 3 reservoir before the 
remainder of the data were taken. This resulted in much smoother 
cooling rates, since the sample chamber would no longer “track" 
each small excursion of the He a reservoir in temperature, but 
rather respond to its general downward trend. This improved 
cooling was reflected in levels of concentration stability that 
were formerly achievable only by the use of heaters. The 
remainder of the data was therefore taken without employing the 
heaters at cooling rates of less than 4 mK/Hr. 

The sensitivity of our results to heat input to the sample 
chamber was a cause of some concern. Since the relatively large 
size of the sample chamber precluded the use of any sort of heat 
shield, we feared that the heat leak to the sample chamber might 
adversely affect our results just as had the use of the heaters, 
in order to investigate this possibility, runs were conducted in 
both warming and cooling modes, presumably affecting both the 
magnitude and the pattern of the heat input to the sample chamber. 
The results were identical to within the accuracy of the experi¬ 
ment. We therefore concluded that our results were not bemg 







influenced by any heat flush effect. 


These results are displayed in Fig. 20. The horizontal axis 
displays the deviation of the concentration from its critical value 
in mole - %. The vertical axis displays the distance below the 
trkritical interface in milimeters. The solid curves are plots 
of eq, 5.1 with the "best fit” values for A, 6, X c and z c . The 
2.5 atm. data and the 9,9 atm. data are shown in Fig. 20(a), 
demonstrating the change in gravitationally induced concentration 
gradient between those two pressures. The 16.7 atm. data are 


shown in Fig. 20(b) for clarity, as it would otherwise overlay the 
2.5 atm. data almost exactly. These results are shown in tabular 


form below: 


Table 5.1 




cm_\ 


6 



pressures. In all cases values of X 
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analysis. Since A is a very strong 
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As can be seen, these 


ie results are consistent with 6-2 at all 
i values of X c and z c that resulted in a 
within the limits determined by graphical 


strong function of 6, we have also 


itude factor when 6-2, as a measure of 
ation gradient. The relatively large 
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Figure 20 

a) Data taken at 2.5 and 9,9 atmospheres plotted o nX-z 
coordinates. Note that the concenti*ation gradient is 
reduced by nearly a factor of two at the higher pressure. 

b) Data taken at 16.7 atmospheres on the same coordinates. 
This graph is shown separately for clarity, as it almost 
exactly overlays the 2.5 atm data. 
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uncertainties at 9.9 atm. are a result of the reduction in the size 
of the effect at that pressure. It is approaching the limit of reso¬ 
lution of our apparatus. No values of X c were determined, since 
only differences in concentration were measured, in addition, 
the exploration of various pressures close to 9.9 atm in an 
attempt to achieve a favorable value of (X -X c ) led to the 
conclusion that the maximum value of X c occurs for a pressure 
of 9.9 ± 0.4 atm. 







VI. Conclusions and Summary 
A. Comparison with Theory and Previous Experiment 

All of our values of 6 are consistent with 6=2, This is in 
agreement both with all the theoretical models considered and 
with the SVP experimental results of Zimmermann 13 et al, and 
Goe liner and Meyer. 13 That the value of 6 seems independent of 
pressure in a system all of whose other parameters vary strong¬ 
ly with pressure is in support of the principle of universality. 

The magnitudes of our measured concentration gradients are 
in agreement with the one other measured value, due to Watts 14 , 
at SVP. Again, assuming 5=2, his data (one point) would lead 
to a value for A z of 1.6±£®. An extrapolation of our values to 
P=0 t although rather uncertain, does not disagree with this 
value. More importantly, correction routines generated from 
our results correctly predict the deviations in concentiation 
observed by Rohde 19 and Johnson 16 in their explorations of the 

behavior of the TCP under pressure. 

A much more stringent test of the magnitudes of the above 
measured concentration gradients results from the calculation 
of (a*/8A) Tp by eq. 3 .18. Determining dX/dz from our values of 
A and 5 by eq. 5.1 and estimating {dv/dX)^ by the methods out¬ 
lined in Appendix D, we have calculated (3A/aX) TP as a function 
of {X-X c ) for the pressures we have employed. These i elation 
ships are plotted in Fig. 21 for values of (*-**> for which we 
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Figure 21 


A comparison with previous experimental work. Our calculated 
values of (BA/8 X) tp compared with those of Goellner and Meyer/' 
The shaded area indicates the uncertainties in the data by Watts 1 * 1 
assuming 6 = 2. 
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« hive data, along with the one resulting Iron, our estimate of 
Walt's" SVP concentration gradient. The data of GoeUner and 
Meyer 15 are displayed as well. As can be seen, our data appears 
entirely reasonable when compared to these previous results. This 
especially true since Watts believes that systematic errors re¬ 
sulted in his underestimation of the concentration gradient near 
the TCP. It should be noted that our data covers another order 
of magnitude closer to the TCP. 

The behavior of the concentration gradient and hence the con¬ 
centration susceptibility under pressure deserves special attention. 
Although our density of points is low, the decrease in magnitude of 
(0X/3A) T1 , between 2 .5 and 9.9 atm. and its subsequent increase as 
the pressure is raised to 16.7 atm. strongly suggests that its 
behavior is analogous to the behavior, as a function of pressure, 
of both the T=0K solubility of He 3 in He 4 and the value of X c . 

More measurements at intermediate pressures will be necessary 
to confirm this. As mentioned in section 3D, we examined the 
Blume Emery and Griffiths model in hopes of shedding some light 
on this behavior. The pressure dependence thus predicted is 
exactly opposite to the one measured. The theoretical values of 
decrease as the pressure is increased from SVP to 
9.9 atm. and then increase up to the solidification point. The 
concentration susceptibility assumes a reasonable value on y 
the if s0 case. Furthermore, to achieve values of X c close 
the 9.9 atm value, both 3Z and (9A/3*) T p become negative. 







A negative value of (»&/BX)„ is indicative of physically un- 
stable solutions. 

B. The Pressure Dependence of the Trieritical Point 

The magnitudes of the concentration gradients measured here, 
as well as the determination of that pressure for which X c has a 
maximum value were used to correct the data of Johnson 18 and 
Rohde 19 near the TCP. The results of this correction are shown 
in Fig. 22, which demonstrates the motion of the TCP under 
pressure on an expanded scale. Where the gravitational gradient 
is large, corrections of .004 in concentration and 3 mK in tem¬ 
perature were common. 

C. Conclusion 

The exploration of the spatial distribution of concentration 
near the TCP is a powerful tool for the determination of both the 
magnitude and critical behavior of tricritical phenomena in He 3 -He 4 
mixtures. A careful determination of a large numboi of points at 
any pressure could determine the critical exponent to within a 
few percent. The spatial resolution also helps to gi'e a much 
clearer picture of the dynamics of interfacial motion. The values 
of the critical exponent and the concentration susceptibility cal¬ 
culated from this data are in good agreement with theory an 
previous experiment. The variation of concentration susceptibility 
with pressure is both marked and remarkable suggesting a 
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Figure 22 

The pressure dependence of the tricritical point. An enlarge 
ment ol the significant portion of Fig. 2. Corrected data by 
Johnson 18 and Rohde 19 , 
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connection between the T = OK solubility of He 3 in He 4 , the value 
ofX c and the value of (3X/&A) TP . it should be noted that the 
van Leeinven and Cohen model, while giving incorrect results 
for the pressure dependence of both the value of and the 
T=0K solubility of He 3 in He 4 , does correctly predict the rela¬ 
tion between them. That is, as one of the quantities increase, 
so does the other. The calculated size of the gravitationally 
induced concentration gradient can be extrapolated to intermediate 
pressures with sufficient accuracy to permit corrections to 
existing data, enabling a more precise determination of behavior 
near the TCP. 









Appendix A 


If we choose X 0. (0, an examination of the molar volume vs 
concentration curves 1 *, Fig. 23, reveals that at SVP the molar 
volume is 33.4 cc and at 17,7 atm the molar volume is 27,0 cc 
This, in turn, gives us a compressibility of 0.012 atm" 1 . 

In a sample chamber 10 cm high, the density of the mixture 
{.UStf/cc) leads to a pressure head of 1 x 10' 3 atm at the bottom 
of the chamber. Thus, the variation in the density of the sample 
over this height is: 

— = 1.2 X 10” s (Al) 

P 

Since the filling of the sample chamber corresponds to a fre¬ 
quency shift of ~ 1.5 x 10 5 Hz, the above variation in density 
corresponds to a change of less than 2 Ilz over the entire height 
of the sample chamber. This is near the limit of resolution of 
the present apparatus. 


*The origin of these curves is 


discussed in Appendix D 













Appendix B 


A number of Maxwell relations can be derived from the 
various thermodynamic potentials appropriate to this system. 

For example, choosing the molar Gibbs potential o£ eq {3 7)’ 

G~u- Ts + Pv 

with: 

dG*-sdT+vdP+AdX 
we can see that: 

©„-* - (§)„-* 

So: 

8 a G fBv \ . ( & Z G \ 

8 Xdp \ &A' anc i BP8X ) 

but, the above second partial derivatives must be equal, giving: 

(to\ _/8A\ (B3) 

Wtt UP/xx 

By choosing other potentials and other derivatives a whole 
series of these relations can bo built up: 


-( 


8A \ 
8P )„ 


(Bl) 


(B2) 
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[»*}„ \»p/ Tlt 

(B6) 

(“1 *(-) 

Win \exl„ 

(B7) 

(!£} -f-') 

\8A/ S , \»»/, 4 

(B8) 

(SL-(g)„ 

(B9) 


(BIO) 

this work, only relation B3 is employed. 







Appendix C 


This is a derivation of the critical behavior predicted by the 
Blume, Emery and Griffiths^ 0 model, solved in the mean, field 
approximation, for general 3 and k parameters. For the super¬ 
fluid case, (A/#0), it is possible to divide eq. (3.55) by eq. (3.56) 
to obtain 


l-X-Mcoth^Af 


Expanding the hyperbolic cotangent and setting <3 = 1: 




( 0 a /) 3 2{m /)$ 

45 + "945 


(Cl) 


(C2) 


In order to go over the reduced variable x -A' c - X, we note 
that Blume, Emery and Griffiths have calculated: 

__ L a _<C3) 

0 2X+Z’ Pc 1 -X c 2JC+1 


So 


X 



(J’Af* 2 g 5 M 6 
45 + 945 


This series can be inverted to give: 


0*Af 2 ss30x 


+ 


r* 2 

15 


+ 0{x 3 ) 


(C4) 


(C5) 


We can re-write eq. (3.55) as: 















exp{j?[A -3C(1-A")]} - ~L- cosh 


(C6> 


K we divide this equation by an identical one in which A and 
X have taken on their critical values, (AT = 0), and take the 
reciprocal, we arrive at: 


- A)} <C7) 

on 

«A C - A) — PX(X C - X) + +ln(l -X) - ln(A') 

- ln[ cosh(3Af] . (C8) 

Substituting (X c -^r) for X in the above equation and expanding 
the last three terms leads to: 


«A c .A)a.p3C* + ln( 55 ^ rT ) + T A^- •* 

J^r 2 

+ ln(l-X c )-ln(X c ) + ^ + |^ + ' * * 

(j3Af) 2 mi .... (C9) 

“ “2~ + 12 


Substituting our value for j3 2 A/‘ i from eq. (C5) into (C9) a 
gathering together like terms, leads, after much cancellation, 
to: 


.r480JC 5 4 1680K 4 4 2368^+lT76^+^M4K+27 j 
A c -A=* 2 --- 120(2K + 1) 3 

(CIO) 


















Thus, fi = 2 and the amplitude of (aA/a*) Tp and hence the con¬ 
centration gradient is a strong function of sc. 

For the normal fluid case, (Jl/ - 0), we can simply look at 
eq. (3.55} with Af - 0. This gives: 

exp{^i} =exp ( 03C(1 -X) J- (Cll) 

Proceeding analogously to the derivation above, we arrive at 

e*p{<l(A c - a)} = eX p{«x-x c )} (C12) 


Taking the logarithms of both sides of this equation and, as 
before, substituting (A' c - at) for X we get: 


which, after some cancellation, leads to: 


■ X r " 2(1 -X c f 


(C13) 


A c -A 


|x + 


(2X - 1)(2K+ 3J 2 
8(23£+l) 


+ • 


(C14) 


Thus, in the normal phase t = 1 and the magnitude of the con 
centration gradient is independent of K to fii st older. 















Appendix D 


The concentration sensitivity of the TDO can be derived as 
follows: the oscillator frequency is a measure of the dielectric 
constant, €, of the material in between the plates of the tank 
circuit capacitor. This frequency is given by: 

1 

^ = VI(cC + C') * D] 


Where L is the total inductance of the tank circuit, C is that 
part of the total capacitance which "sees' 1 the sample and C' is 
the remaining "stray” capacitance. A rearrangement of terms 
gives: 


' = 4* S LC/‘'(I") 


(D2) 


Our first task is to determine the magnitude of C'. From (Dl) 
it can be seen that if the sample chamber is empty, the frequency, 
/t» will be given by: 


_ 1 

^ 0= 2vVLB€ 


* C' 
where B = 1 + -rr 


(D3) 


Substituting (D3) into (D2) results in 


c-1 



(D4) 


In order to evaluate B the sample chamber was first cvacu 
■ad then filled with He*. The frequencies before and after fillmg 












were noted. The dielectric constant of pure He 4 at SW is known 
Irom data by Kerr and Sherman 37 and thus B may be cal¬ 
culated. The results of this procedure were: 

/ 0 = 11,544,096 ±5 

/ = 11,395,760 ± 5 

€ = 1.0573 ± .00005 

B=2.19± .01 

Once B is known, changes in frequency can be equated with 
changes in dielectric constant. The dielectric constant can be 
related to the molar volume of the sample by the Clausius- 
Mosotti relation: 

Vm = ll£i±i <D5> 

3 €-1 

where P is the molar polarizability of helium, which is, for oui 
purposes, the same for He 3 and He 4 ,* 7 and Vm is the molar volume 
of the sample. Finally, the molar volume of the sample can be 
related to its concentration by generating a set of Vm vs. A' 
curves at various pressures. Such a set of curves is sho 
Fig. 23. The data for pure He 4 is by Kerr and Sherman 37 and 
Keesom 58 and the data for pure He 3 is by Sherman and Edeslcu y 
*>ataat intermediate points at SVP is due to Kerr 40 . The curves 
at higher pressures were arrived at from analysis of th 











Figure 23 


Molar volume vs. concentration curves for He 3 -He" mixtures. 
Straight lines indicate ideal solution behavior. SVP data is due to 
Kerr. 40 The data Cor X = 0 is by Kerr and Sherman 37 and Keeaom . 
The data for X = 1 is by Sherman and Edeskuty 39 . Intermediate 
points at elevated pressure were arrived at by analysis of the 
data of Johnson 18 and Rohde 19 . 
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Johnson 1 ' and Rohde* \ who measured the molar volumes at two 
intermediate concentrations exactly , and extrapolated the re¬ 
mainder of the curves on the basis of self-consistency arguments. 
We can, thus, in principle, determine the absolute concentration 
of our sample by the method outlined above. 

In fact, during the course of this experiment, we were in¬ 
terested only in changes in concentration relative to the average, 
and so a simpler procedure is possible. By differentiating (D5) 
we can obtain: 


dVrn ZVnt 


(D6) 


rfc (c-1)(£ + 2) 


Then a combination of (D4) and (D6) yields: 


dVm 

rfc 



3 Km 


(D7) 


but from (D4) 


df 


2 B/g 
/ 3 


(D8) 


So: 


9 ^ 

Since the molar volume of He 4 at SVP is 27.52 cc we can 



evaluate {D9) to give 
















■=1,83 X 10‘ 4 at that pressure . 

df 


(DIO) 


Then the concentration sensitivity of the oscillator may be 
found thusiy: 

dX ^ dx dVm 

df = dVm df (D11) 

Since over the course of a run the changes in X and/are very 
small compared to their values — no more than 1%—we can 
assume that; 




dX 

dVm 


dVtn 

~df 


4f 


(D12) 


Values of dX/dVm can be estimated at any pressure from the 
values of Vm for pure He 3 and pure He' 1 at that pressure. As can 
be seen from the curves of Fig. 23 the approximation that. 


dX 1.0-0.0 

dVm = Vm(He 3 ) - Vm(He 4 ) 

is quite good in the range of interest - 0.67 ^X - 0.70 as the 
tangent and chord of the curves are nearly parallel in that range 
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